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We show that the recently proposed S = 2 bond model 
with orbital degrees of freedom for insulating V2O3 not only 
explains the anomalous magnetic ordering, but also other 
mysteries of the magnetic phase transition. The model con- 
tains an additional orbital degree of freedom that exhibits a 
zero temperature quantum phase transition in the Ising uni- 
versality class. 

PACS numbers: 71.30. +h, 75.10.-b, 75.50.Ee 



The metal-insulator transition in vanadium oxide 
(V2O3) has long beenJiailed as a physical realization of 
the Mott transitionErEl. However, the insulating phase 
displays essential deviations from the standard Mott in- 
sulator with only spin degrees of freedom. For example, 
while it undergoes a transition from paramagnetic to an- 
tiferromagnetic (AF) insulator with decreasing temper- 
ature, the magnetic ordering pattern observedaTJ (hith- 
erto called RS) is not that of the simplest two-sublatticc 
Neel state (hitherto referred to as AS). Some time ago, 
Castellani, Natoli, and Ranninger (CNR)Q proposed a 
model in which the low energy degrees of freedom for each 
V ion are a spin (5 = 1/2) coupled to doubly degener- 
ate orbital. They showed that within certain parameter 
regime, the ground state exhibits the RS magnetic or- 
dering pattern together with orbital ordering. Recently, 
resonant x-ray scattering experiments!! have observed or- 
bital ordering in V2O3 with an ordering wavevector con- 
sistent with one of the ordered state (so-called RO state) 
obtained by CNR. 

The CNR model relies on one out of the two electrons 
on a V ion forming a spin singlet bond with its coun- 
terpart on the neighboring V ion along the c-axis, leav- 
ing one electron occupying doubly degenerate orbitals on 
each V ion. This has recently been criticized for not re- 
specting the strong on-site coulomb repulsion and Hund's 
rule couplingpt2l. This criticism is supported by polar- 
ized soft x-ray experimentsQ, and strongly suggests that 
the effective spin on a V ion should have 5 = 1 instead 
of 5 = 1/2. Recently, a 5 = 1 model incorporating the 
orbital degree of freedom has been introduced!!!. In this 
model V-V pair along c-axis is locked ferromagnetically 
into a total spin 2 "bond" with a double orbital degen- 
eracy. These bonds occupy a lattice that is topologically 
equivalent to a cubic lattice (Fig. 1). Virtual in-plane 
hopping in the honeycomb planes couple such bonds to 
their neighbors. For Hund's rule coupling and hopping 
constants consistent with reported/expected values, it 



was shown that the ground state has RS magnetic and 
ferro-orbital (FO) order. The FO ordering was shown to 
be also consistent with the resonant x-ray experiments. 
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FIG. 1. Mappipgi-ef bonds on the corundum lattice to sites 
on a cubic latticetillii The positions of V ions are shown by 
squares and the centers of the bonds by circles. The arrows 
represent the spin of V-ions in V-V pair in the RS state. 

In addition to the anomalous RS magnetic order, there 
are other mysteries associated with the V2O3 insulator!! 2 ]: 
I) the magnetic and orbital ordering occur at the same 
transition temperaturea, in contrast to the Jahn- Teller 
transitions observed in other transition metal-oxides like 
MnO; II) the transition is strongly first orded!3; and III) 
above the transition temperature, neutron scattering ob- 
serves a relatively broad peak in the susceptibility for 
wavevector of AS ordering but no features for the RS 
ordering wavevectonij. 

In thus letter, we will show that the 5 = 2 bond 
modelEU gives satisfactory explanation for the three mys- 
teries above. In addition, we will elucidate on a new or- 
bital degree of freedom, which can lead to a second orbital 
transition at a lower temperature. The critical tempera- 
ture of this second orbital transition goes to zero as the 
intraplane hopping approaches some critical value, giv- 
ing rise to a zero temperature quantum phase transition. 
This second transition is well described by the transverse 
field Ising model. 

We begin by briefly recalling the steps leading to the 
effective 5 = 2 bond Hamiltonian. The large on-site, 
same-orbital coulomb repulsion (U « 5eV) and Hund's 
coupling (J w 0.6 — leV)&E3 constrains the Hilbert space 
of each V ion to that of 2 electrons with total spin 5=1, 
hence singly occupying two of the lower three t2 g d- 
orbitals favored by crystal field. The remaining 9-fold 
degeneracy of a single ion is then lifted due to virtual 
hopping. Let tf • be the various hopping integrals between 



1 



nearest neighbor (nn) V ions, where S = a,b,c,d stands 
for the direction of the bond (Fig. 1) and = 1,2,3 
for the orbitals e g %, e g 2, ai g respectively. Since £33 is by 
far the largest, the strongest effect is due to virtual hop- 
ping between vertical pairs of V ions. When combined 
with J large, this further reduces the Hilbert space of a 
single pair to that of total spin S — 2, with an occupation 
of (e g i, e g 2) by the two electrons on one V ion and (e g , 
dig) on the other. The orbital configuration proposed 
here is qualitatively consistent with the polarized x-ray 
absorption spectra of V2O3EI. The basis for such a V-ion 
pair can then be written as \a z Tfj). Here a z — —2, —1, ...2 
is the spin state, r = +1 (—1) if e g \ (e g 2 ) is occupied 
on the ion with (e ff , a% g ) occupation, and fi = +1 (— 1) 
denotes whether the top V ion has the (e 9 i, e g 2) or (e ff , 
dig) occupation. 

Within this restricted Hilbert space, we can now obtain 
the effective Hamiltonian by considering the lifting of the 
degeneracies due to the effects of the other tfj's. The ef- 
fective Hamiltonian is expressed by defining on each ver- 
tical pair a spin 2 operator er, a psuedospin 1/2 operator 
r, and yet another psuedospin operator /it, which act on 
the a z , r, and fi degrees of freedom respectively. The ef- 
fective Hamiltonian must respect global SU(2) symmetry 
in spin space, Z3 symmetry in r space, and Z2 symmetry 
in fj, space. The Z3 Potts symmetry in r reflects the C3 
rotational symmetry of the corundum lattice about the 
c-axis while the Ising Z2 symmetry in fj, reflects the global 
"top-bottom" inversion symmetry of the V-ion pairs. We 
consider only processes where after two (virtual) hops, 
the system is restored into the restricted Hilbert space. 
With no loss of generality, we assume for i?- that only 
tn and £23 are non-negligible. 

At temperature « 155K, x-ray experiment haSpShown 
that all V ions remain equivalent in the AF phasecJ, im- 
plying no symmetry breaking of the Z2 /i symmetry at 
the magnetic transition. If we further neglect fluctuation 
effects, effectively setting \i z = in the effective Hamil- 
tonian, fi becomes decoupled from the a and r. deg rees 
of freedom, and the Hamiltonian is simplified totirE-3 

-ff = a ^ S, ; • Sj +ai 2J • Sj(ri ■ Sy- + tj ■ Ufj) 
+a 2 Tj • UijTj ■ + b 2 Si ■ Sj (n ■ UijTj ■ Uij) 

(ij) (ij) 

where Sy = Si, S2, S3 for i,j nn in the x, y, z directions 
respectively. The Sj's are unit vectors in the x — z plane 
of the r space, with S3 = z, while Si and S2 are ro- 
tated from S3 by 120° and 240°. The coupling constants 
ao etc. depend on J and 7 = (ti 1/^23) ■ The model 
contains three types of local order parameters: a vector 
order parameter M, =< Si >, a Pott's type order pa- 
rameter rf =< Ti ■ n a >, and a tensor order parameter 
Qf =< SiTi-n a >. Assuming up to four sublattice order- 
ing in these three order parameters on the equivalent cu- 
bic lattice with wavevectors k x — k y = 0, n, k z = 0, n, we 



approximate H by a single-site Hamiltonian using usual 
mean field theory (MFT) to find the values of Mj, rf, 
and Qf that minimizes the free energy. The phase dia- 
gram for zero temperature thus obtained is shown in Fig. 
2. 
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FIG. 2. T — phase diagram of the S — 2 bond model. 
U' = U - 2J. Phases are 1) FORS, 2) AOFS, 3) FORS', 
4) AOAS, 5) AO z RS, and 6) AO^RS' (subscripts after AO 
indicate direction of AO ordering in r space). Dashed line 
shows where ao = 0. Dotted line shows the boundary between 
first order (below the line) and second order transitions for 
the FORS phase. 

The phase diagram can be understood qualitatively as 
followsEJ'Ell The key dependences of the various cou- 
pling constants on J and 7 are that ao changes from 
> to < as J is increased, ai < everywhere, and 
0-2^2 = for 7 = and > otherwise, with 02 > 4J>2. 
The ao term favors AS or ferromagnetic spin (FS) order- 
ing while the 02 and 62 terms combined favor antiferro- 
orbital (AO) ordering. The ai term in H appears as an 
explicit symmetry breaking field on r. However, for the 
AS or FS states or any other state whose spin correlations 
< Si ■ Sj > do not break the lattice C3 symmetry, this 
field equals zero when summed over all bonds, and the 
ai term gives no contribution to the energy. Thus, close 
to the boundary where ao changes sign, the spins should 
order magnetically in a way to break the C3 symmetry 
to take advantage of the ai term. This accounts for the 
states RS and RS which differ from each other by having 
1 ferromagnetic bond and 2 antiferromagnetic bonds on 
the honeycomb plane for the RS and vice versa for RS . 
Choosing the ferromagnetic bond for RS and antiferro- 
magnetic bond for RS to be along the a direction of the 
honeycomb planes, these magnetic orderings set up an 
uniform field on r in the +z and — z directions respec- 
tively, thus for small 7 causing (tj Z ) to become uniformly 
non-zero (FO phase). 

Of these phases, the one that is consistent with the 
magnetic and orbital ordering observed by neutron scat- 
tering and resonant x-ray scattering experiments is the 
FORS phase. Below we shall focus on this phase and 
study the finite temperature phase transition of the 
model to explain the three mysteries of insulating V2O3. 

I) In V2O3, the magnetic transition is accompanied by 
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orbital ordering. In our MFT, we indeed find the spin 
and orbital to order at the same temperature. This re- 
sult is robust, and can be explained more generally as 
a consequence of the form of the Hamiltonian together 
with FORS type ordering. Since the RS ordering pro- 
duces a symmetry breaking field on r, any RS ordering 
must necessarily be accompanied by non-zero (rj)'s. The 
converse of this is not true, and orbitals can in principle 
order without spin ordering. However, the only term in 
H that can cause FO ordering is the a\ term which van- 
ishes without RS ordering!!! Thus, we conclude that for 
FORS (or FORS ) ordering, magnetic and orbital order- 
ing must occur together. This is not true for the other 
ground states. For AOAS or AOFS, as T is lowered, mag- 
netic ordering can first occur without orbital ordering or 
vice versa. For AORS and AORS , magnetic ordering 
must be accompanied or preceded by orbital ordering, 
but orbital ordering can occur without magnetic order- 
ing. 

We should note that the observed monoclinic distor- 
tion below T c plays no role in our theory. Indeed, we 
believe it to be a by-product of the FORS transition 
breaking the C3 rotation symmetry of the honeycomb 
planes. 

II) The magnetic transition is strongly first order, with 
an entropy jump of w ks In 2 per V ion. The size of en- 
tropy jump suggests that the first order transition should 
not be induced by critical fluctuations, but should be 
mean field in origin. Indeed, our MFT shows a first 
order transition for the paramagnetic-FORS transition 
for small 7 (Fig. 2). It should be remarked that while 
the symmetry of H allows third order invariants in the 
Landau-Ginzburg free energy functional of the form r 
and MQr, they are not the cause of the first order tran- 
sition. This is because the FO susceptibility at T c is small 
since the FO ordering is not favored by the bare orbital- 
orbital coupling of H . Indeed, for values of 7 such that 
the transition into FORS phase is continuous, r ~ MQ 
below but close to the transition. Instead, the first or- 
der transition is due to the presence of a negative M 2 Q 2 
quartic term. 

While MFT gives a first order transition, the calcu- 
lated entropy jump and maximum value of 7 are rather 
small, especially when compared to the actual experi- 
mental value of w fc B In 2 per V ion. We believe this to 
be a result of ignoring nn AF correlation in the param- 
agnetic state in a MFT for the FORS transition. This 
correlation is seen to be sizeable experimentally (see III 
below), and serves to stabilize the paramagnetic state. 
We estimate the nn correlation using high temperature 
expansion to second order in ao and find the entropy 
jump and maximum 7 to be significantly increased with- 
out much reduction in T c . A schematic of the free energy 
functional without and with nn AF correlations included 
is shown in Fig. 3. 

III) Given the magnetic RS ordering below T C) one 
might expect peaks at the three symmetry related 
RS wavevectors above T c , and that these peaks would 




FIG. 3. Solid line shows schematically the MF free energy 
as a function of the FORS order parameter. Dotted line shows 
the free energy when nearest neighbor AF correlations are in- 
cluded for the paramagnetic state. Inclusion of such corre- 
lations stabilizes the PI state and lowers the transition tem- 
perature slightly from to T c while significantly increasing 
the order parameter and entropy jump. 

sharpen as theptemperature is lowered towards T c . Neu- 
tron scattering^ above T c , however, shows a broad peak 
in AS wavevector but no feature in the RS wavevectors. 
Moreover, as T decreases, the broad peak increases in in- 
tensity but does not sharpen appreciably, indicating that 
nn spins are increasingly AF correlated but the correla- 
tion length remains of order of nn. This seemingly sur- 
prising feature can be explained qualitatively as follows 
in our theory. Above T c , there are fluctuations in both 
the AS and RS states. The nn spin correlation of AS is 
of course AF. While the correlation of RS for a particular 
symmetry broken state has a pattern with ferromagnetic 
in one nn bond and AF in other two nn bonds in the hon- 
eycomb plane, the nn correlation when averaged over the 
three RS states also is AF. For long distance correlations, 
however, the signs of correlations due to fluctuations of 
AS and RS states can be competing. As a result, these 
spin correlations remain weak. 

FO ordering exists only for RS or RS ordering. Also, 
the simultaneity of magnetic and orbital ordering, and 
the first order phase transition depend crucially on FORS 
ordering. Thus, the phase transition phenomenology can 
be changed by changing the type of magnetic ordering. 
This can be achieved by applying an uniform magnetic 
field which with increasing strength will change the pre- 
ferred magnetic ordering from RS to FS. This will in 
turn change the orbital ordering and the phase transition 
properties. This can be observed by performing resonant 
X-ray scattering in the presence of a magnetic field. 

Let us now return to the y, degree of freedom. At tem- 
perature significantly lower than T c , the FORS ordering 
will be approximately saturated and frozen out, and fi 
becomes the only relevant degree of freedom. The effec- 
tive Hamiltonian thus becomes 

« (y> 
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where A = 2*J 1 tf 3 /(Z7 - 3J), Jy 



cally aligned nn, and J^ 
Ji and J2 are given by 



Ji for ferromagneti- 
J2 for antiferromagnetic ones. 



Ji = 



2:i 



■h 



3J' 

2 +i *—YH-~i) 

(U-3J) U U + 2JJ 24 v 8 ,; 



We recognize this Hamiltonian as the transverse field 
Ising model (TFIM). At T — 0, this model has a quantum 
phase transition from a disordered state for large A to an 
ordered one for small A. The ordering pattern depends 
on the signs of Ji and J2. We see that Ji > always, 
while J2 changes sign as 7 is decreased. For J2 < 0, /1 or- 
dering does not lead to further reduction of translational 
invariance from the FORS state, and no new wavevector 
q will be observed in resonant x-ray diffraction, only in- 
tensity change. For J2 > 0, there will be further change 
in Bravais lattice symmetry, and peaks at new g's will 
be observed. The value of 7 necessary for FORS order- 
ing and in particular to have a first order transition for 
estimated value of J suggests that J2 > 0. 

For the 3D TFIM, its finite temperature transition is in 
the same universality class as the 3D Ising model and its 
T=0 transition in that of the 4D Ising model respectively. 
However, if J2 is close to zero, then the system behaves 
like a quasi-l-D one. Therefore, as J2 varies, it may be 
possible to observe a dimensional crossover. If J 2 is not 
too small, the critical value of A for LRO may be ob- 
tained by MFT as 2Ji +4 | J 2 | = A/ tanh(A/fc B T). This 
implies the critical temperature rises sharply from 
as A is reduced from its T — critical value. The MF 
value of T M should be accurate provided A is neither too 
close to A c nor too small. For A close to A c , Tp is gov- 
erned by quantum critical phenomena, and scales as the 
inverse correlation time, so that up to logarithmic cor- 
rections T M ~ £7" ~ (A c — A). For small enough A, 
T M can become comparable to T c , and the assumption of 
saturated FORS ordering will not be valid. 
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